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Abstract. The Large Scale Structure in the galaxy distribution is investigated using The First Data Release of 
the Sloan Digital Sky Survey. Using the Minimal Spanning Tree technique we have extracted sets of filaments, 
of wall-like structures, of galaxy groups, and of rich clusters from this unique sample. The physical properties 
of these structures were then measured and compared with the statistical expectations based on the Zel'dovich' 
theory. 

The measured characteristics of galaxy walls were found to be consistent with those for a spatially flat ACDM 
cosmological model with Q, m w 0.3 and Qa ~ 0.7, and for Gaussian initial perturbations with a Harrison - 
Zel'dovich power spectrum. Furthermore, we found that the mass functions of groups and of unrelaxed structure 
elements generally fit well with the expectations from Zel'dovich' theory. We also note that both groups and rich 
clusters tend to prefer the environments of walls, which tend to be of higher density, rather than the environments 
of filaments, which tend to be of lower density. 
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1. Introduction 

With the advent of the Durham/UKST Galaxy Redshift 
survey (DURS, iRatcliffe et all Il99flh and the Las 
Campanas Redshift Survey (LCRS, IShectman et alJ 
the galaxy distribution on scales up to ^300 
hr 1 Mpc could be studied. Now these investigations can 
be extended using the public data sets from The First 
Data Release of the Sloan Digital Sky Survey (SDSS DR1, 
IStoughton et alJl2002t fAbazaiian et al.ll2003j) . which con- 
tains redshifts for w 100000 galaxies in four slices for dis- 
tances D < 600ft- x Mpc. 

The analysis of the spatial galaxy distribution in the 
DURS and the LCRS has revealed that the Large Scale 
Structure (LSS) is composed of walls and filaments, that 
galaxies are divided roughly equally into each of these 
two populations (with few or no truly isolated galax- 
ies), and that richer walls are linked to the joint ran- 
dom network of the cosmic web by systems of filaments 
l|Doroshkevich et alJl200oll200l|) . Furthermore, these find- 
ings are consistent with results obtained for simulations of 
dark matter (DM) distributions (see e.g. lCole et al.lll99a: 
Ijenkins et al.ll998(l an d for mock galaxy catalogues based 
upon DM simulations jCole et al.lll99S|h 

The quantitative statistical description of the LSS is 
in itself an important problem. Beyond that, though, the 



analysis of rich catalogues can also provide estimates for 
certain cosmological parameters and for the characteris- 
tics of the initial power spectrum of perturbations. To do 
so, some theoretical models of structure formation can be 
used. 

The close connection between t he LSS and Zel'dovich' 
panca kes has been discus sed by iThompson fc Gregory! 
Hizj) and by lOorti l)l983|) . Now this connection is ver- 
ified by the comparison of the statistical characteris- 
tics of obse rved and simulated walls with theoretic al ex- 
pectations ijDemianski fe Doroshkevicblll999l I2002L here- 
after DD99 and DD02) based on the Zel'dovich t heory 
of nonlinear gravitational i nstability l|ZerDovichl Il97fi 
IShandarin fc Zeldovicfll989lb This approach connects the 
characteristics of the LSS with the main parameters of 
the underlying cosmological scenario and the initial power 
spectrum, and it permits the estimation of some of these 
parameters using the measured properties of walls. It was 
examined with the s imulated DM distribution (DD99; 
IDemianski et al l2000lh and was found that, for sufficiently 
representative samples of walls, a precision of better than 
20% can be reached. 

Effective methods of the statistical description of the 
LSS based on the Minima l Spanning Tree ( MST) tech- 
nique were developed by IDemianski et al.l l|200fj|) and 
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Fig. 1. Four regions of the DR1 sample on the sky. 



iDoroshkevich et all (|2000l l200l[) . who applied them to 
DM simulations and to the PURS and t he L CRS. 
These methods introdu ced in IB arrow et al.1 l|l985j) and 
Ivan de Weveaer 3 lll99lh generalize the popular "friends- 
of-friends" approach. In this paper we apply the same 
approach to the SDSS DR1, a sample from which we 
can obtain more representative and more precise mea- 
sures of the properties of the LSS and the initial power 
spectrum of perturbations. Alternative methods based 
on Minkowski functiona ls a re proposed by, fo r example, 
ISchmalzing et alJ l|l999j) . or lSheth et all (|2002|) 

With the MST technique, we can quantitatively de- 
scribe the sample under investigation and divide the sam- 
ple into physically motivated subsamples of clusters with 
various threshold overdensities bounding them. Then we 
can extract the LSS elements and characterize their mor- 
phology. In particular, this technique allows us to discrimi- 
nate between filamentary and wall-like structure elements 
located presumably within low and high density regions 
and to estimate their parameters for the different thresh- 
old overdensities. The same technique allows us to extract 
sets of high density groups of galaxies and to measure 
some of their properties. 

The analysis of wall-like condensations is most in- 
formative. Comparison of the observed characteristics of 
walls with theoretical expectations (DD99; DD02) demon- 
strates that the observed galaxy distribution is consistent 
with Gaussianity initial perturbations and that the walls 
are the recently formed, partly relaxed ZePdovich' pan- 
cakes. The mean basic characteristics of the walls are 
consistent with those theoretically expected for the ini- 
tial power spectrum measur ed by the CMB obse rvations 
summarized, for example, in ISnergel et all (|2003h . 

In this paper we also analyse the mass functions of 
structure elements selected for a variety of boundary 
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Fig. 2. The radial galaxies distributions in the four sam- 
ples of the SDSS DR1. The selection function Q is plotted 
by solid lines. 



threshold overdensities. We show that these functions are 
quite similar to the expectations of Zel'dovich' theory, 
which generalizes the Press - Shechter formalism for any 
structure elements. In addition, the theory indicates that 
the interaction of large and small scale perturbations can 
be important for the formation of the observed LSS mass 
functions. Our analysis demonstrates that this interaction 
is actually seen in the influence of environment on the 
characterist ics of groups of galaxies. This problem was also 
discussed in lEinasto et alJ l|2003alb|) . 

This paper is organized as follows: In Sees. 2 we de- 
scribe the sample of galaxies which we extracted from the 
SDSS DR1 and the method we have employed to correct 
for radial selection effects. In Sec. 3 we establish the gen- 
eral characteristics of the LSS. More detailed descriptions 
of the filamentary network and walls can be found in Sees. 
4 and 5, respectively. In Sees. 6 and 7 we discuss the prob- 
able selected clusters of galaxies and the mass function of 
structure elements. We conclude with Sec. 8 where a sum- 
mary and a short discussion of main results are presented. 



2. The First Data Release of the Sloan Digital 
Sky Survey 

We use as our obser vational sample the SDSS DR1 
l|Abazaiian et al.l l2003h , which is the first public release 
of da t a from the SPSS l|Fukugita et alJll99fl: lOunn et alJ 
Il998t lYork et aljElOol) . The imaging data for the SPSS 
DR1 encompasses 2099 sq deg of sky. The DR1 also con- 
tains 186,240 follow up spectra, which are available over 
1360 sq deg of the imaging data area. Galaxies are sit- 
uated within two north fields, N1&N2, and two south 
fields, Si &S*2. These regions are plotted in Fig.fJ. 
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Fig. 3. The wedge diagram of galaxy distributions in the 
samples Si is real (top panel) and modified (bottom panel) 
radial coordinates. 



We obtained our SDSS DR1 sample via the SDSS DR1 
Spectro query server 1 . This is a web interface to the SDSS 
Catalog Archive Server. We selected all objects identified 
as galaxies with a 'redshift confidence minimal level' of 
95% and no maximal level. No other constraints in the 
selection were made at this level. 

Our method for detecting LSS depends on having 
largely contiguous regions. Hence we have removed some 
regions and artifacts from the DR1 sample. The follow- 
ing are the RA and DEC areas we masked out from the 
original DR1 query before our analysis: 

• 174 ,l <RA<179' i , -4.0°<DEC<-1.22° 

• 159 ,l <RA<163' i , 1.1°<DEC<4.0° 

• 10 h <RA<50 h , 10°<DEC<20° 

• 300 /l <RA<355 /l , -12°<DEC<-4° 

• 250' 1 <RA<270'\ 52°<DEC<67° 



2.1. Correction for radial selection effects 

In Figure [21 we plot the radial distributions of galaxies 
in all four samples. In Figure we plot the wedge di- 
agramm of observed galaxy distribution for the sample 
Si- As is clearly seen from these Figures, at distances 
D > 400/i~ 1 Mpc this fraction of observed galaxies is 
strongly suppressed because of the radial selection ef- 
fect. Note that this suppression is quite successfully fit 
by curves describing a selection function of the form 

f ga l(D) oc D 2 exp[-0D/i? sei ) 3/2 ], R sel a lQOfc^Mpc (1) 
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Fig. 4. The normalized mean galaxy density in the four 
modified samples of SDSS DR1. 

where D is a galaxy's radial distance and R se i is the selec- 
tion scale l)Baueh fc Efstathioulll993j) . These fits are also 
plotted in Fig. El 

In some applications, like when we want to correct a 
measure of the observed density to a measure of the true 
density, we would like to use equation to correct for the 
radial selection effects after the fact. An example of such 
a case is calculating a group's or cluster's true richness 
based upon the observed number of galaxies it contains 
(Sec. 6 & 7). modified 

In other applications, however, like in searching for 
groups or clusters in a magnitude- limited sample, we want 
to make a preemptive correction for the radial selection ef- 
fects. For example, in a standard friends-of-f riends perco- 
lation algorithm fe.g. lHuchra fc Gellerl ll982). this is done 
by adjusting the linking length as a function of radial dis- 
tance. Here, instead, we employ the rather novel approach 
of adjusting the r adial distances them s elves a s introduced 
and discussed in iDoroshkevich et alJ |2001). Hence, in- 
stead of the measured radial distance, we use a modified 
radial distance, D m d, where 

(D/R sel ) 3 / 2 ]eM~(D/R se i) 3/2 ])(2) 



D 3 md = 2R 3 sel (l-[l 



1 http://www.sdss.org/drl 



The radial variations of the normalized number density of 
galaxies for all samples from Figure|3are plotted in Figure 
As is seen from this figure, the modified radial distances 
for the galaxies suppresses the very large-scale trends. On 
the other hand the relative position of galaxies remain 
unchanged and the smaller scale random variations in the 
density are emphasized. The wedge diagramm of modified 
galaxy distribution for the sample Si is also plotted in 
Fig.0 

This correction does not change distances at D < R se i 
and is more important for the more distant regions of our 
samples (D > 350/i _1 Mpc), which contain only ~20% of 
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all galaxies. Thus, in the following analyses, we apply this 
correction only to the separation of the high and low den- 
sity regions in the deeper samples. Of course, it cannot 
restore the lost information about the galaxy distribution 
in these regions, but it does help compensate for the strong- 
drop in the observed galaxy density at these distances. It 
also allows one to apply the standard methods of investi- 
gation for the full catalogues with a depth of 600/i _1 Mpc. 

3. General characteristics of observed large scale 
structure 

To characterize the general properties of the large scale 
spatial galaxy distribution we use the Minimal Spanning 
Tree (MST) technique applied to both directly observed 
samples of galaxies and to samples corrected for the selec- 
tion effect. 

In our analysis here, we consider the four fields, plotted 
in Fig.rU at the distance D < 420/i~ 1 Mpc with 



N gal = 79 183, 



Igalj 



10 _2 Mpc" 



(3) 



where N ga i and {n ga i) are the total number of galaxies 
and the mean density of the samples. This sample contains 
«85% of all galaxies with moderate impact from selection 
effects. The numbers of galaxies in the separate fields are 

- N\, the northern sample (35 520 galaxies) 

- Ni, the northern sample (21983 galaxies) 

- the southern sample (11225 galaxies) 

- 5*2, the southern sample (10455 galaxies) 

3.1. The MST technique 

The M ST tec hnique was first discuss ed bv iBarrow et alJ 
l|l985|) and bv lvan de Wevgaertl l|l99ll) . The MST is a con- 
struct fro m graph theor y, originally introduced bv lKruskaJ 
l)l956l) and lPri inl lll957h . which has been widely applied in 
telecommunications and similar fields. It is a unique net- 
work associated with a given point sample and connects all 
points of the sample to a tree in a special and unique man- 
ner which minimizes the full length of the tree. Further 
definitions, exa mples, and applicat ions of this approach 
are dis cussed in lBarrow et alJ l|l985|) and lvan de Wevgaertl 
More reference s to the m athematical results can 
also be found in Ivan de Wevgaertl l)l991^ . 

One of earliest uses of MST approach in the s t udy o f 
large-scale structure was that of iBhavsar fc Lind <|l988|h 
who successfully applied it to extract filamentary struc- 
tures from the original CfA Redshift Survey. Its ap- 
plications for the quantitative description of observed 
and simulated catalog u es of galaxies were discussed i n 
iDemianski et alJ l|200fl) : iDoroshkevich et alJ |200o[l200ll) . 

One of the most important features of the MST tech- 
nique is generalization of the widely used "friends-of- 
friends" approach. It allows one to separate all clusters, 
LSS elements, with a given linking length. In spite of 
the very complex shape of the clusters, the linking length 



defines for each two points the local overdensity bound- 
ing the clusters with a relation familiar from "friends-of- 
friends" algorithms l|Huchra fc Gellerlll982h : 



Sthr = 3/[47r(n ffai )rf nfc ] 



(4) 



Further on it allows one to obtain characteristics of each 
cluster and of the sample of clusters forming the LSS with 
a given overdensity. Further discrimination can be per- 
formed for a given threshold richness of individual ele- 
ments. 

Here we will restrict our investigation to our results 
for the probability distribution function of the MST edge 
lengths Wmst(1) and to the morphological description of 
individual clusters. The potential of the MST approach is 
not, however, exhausted by these applications. 

3.2. Wall-like and filamentary structure elements 

With the MST technique we can demonstrate that the ma- 
jority of galaxies are concentrated within wall-like struc- 
tures and filaments which connect walls to the joint ran- 
dom network of the cosmic web. The internal structure 
of both walls and filaments is complex. Thus, wall-like 
structures incorporate some fraction of filaments and both 
walls and filaments incorporate high density galaxy groups 
and clouds. In particular, clusters of galaxies are usually 
situated within richer walls while groups of galaxies are 
embedded within filaments. In spite of this, the galaxy 
distribution can be described as a set of one, two and 
three dimensional Poisson-like distributions. Naturally, a 
one dimensional distribution is more typical for filaments 
while two and three dimensional ones are typical for walls 
and groups of g a laxies , respectively. As was shown i n 
Ivan de Wevgaer 3 l|l99lt) iBurvak fc Doroshkevich! l|l99fi(L 
a Poissonian distribution of galaxies within the LSS ele- 
ments successfully reproduces the observed 3D correlation 
function of galaxies. 

These result show that the probability distribution 
function of MST edge lengths (PDF MST), WmstQ), 
characterizes the geometry of the galaxy distribution. For 
ID and 2D Poissonian distributions, typical for filaments 
and walls, Wmst{1) is described by the following expo- 
nential and Rayleigh functions, 



W MS t(1) = W e (l) = (ly 1 ex P H/(0) , 
W M st(1) = WrU) = 2l/(l 2 ) cxpH7</ 2 )) 



(5) 



These PDFs remain valid for any ID and 2D distributions 
when the galaxy separation is small as compared with the 
curvature of the lines and surfaces. Comparison of mea- 
sured and expected PDFs MST allows one to demonstrate 
the existence of these two types of structure elements 
and to make approximate estimates of their richness. Let 
us recall, however, that for the galaxy groups embedded 
within filaments, 2D and 3D Poissonian distributions are 
observed and, so, in this case we cannot see a purely ID 
distribution. For walls this effect is less important because 
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Fig. 5. PDFs of MST edge lengths in redshift space aver- 
aged over four samples are plotted for the full sample (top 
panel), HDRs (middle panel) and LDRs (bottom panel). 
Rayleigh and exponential fits are plotted by thin solid and 
dashed lines. 



filament s were proposed and tested in our previous publi- 
catio ns l|Demiahski et al.ll200(i iDoroshkevich et alJl200(]| 
200l|. The first step is to make a rough discrimination be- 
tween the high and low density regions (HDRs and LDRs). 

Such discrimination can be easily performed for a 
given overdensity contour bounding the clusters and a 
given threshold richness o f individual elements. Following 
IDoroshkevich et "all l)200lh . in all four samples with D < 
420/i~ 1 Mpc, wall-like high density regions (HDRs) were 
identified with clusters found for a threshold richness 
Nthr = 40 and a threshold overdensity contour bound- 
ing the cluster equal to the mean density, 8 t hr = 1- These 
samples, Ni,N 2 ,Si,S 2 , of HDRs contain 49%, 47%, 51% 
and 47% of all galaxies. The samples of low density regions 
(LDRs) , which are occupied mainly by filaments and poor 
groups of galaxies, are complementary to the HDRs in 
that the LDRs are simply the leftovers from the original 
total samples after the HDRs have been removed. 

In Figure^ (middle panel) the Wmst(1) plotted for the 
HDRs is very similar to a Rayleigh function, thus confirm- 
ing with this criterion the sheet-like nature of the observed 
galaxy distribution within the HDRs. As before, the error 
bars show the scatter of measurements for the four sub- 
samples. For 90% of objects we have 



it only distorts the 2D distribution typical for such LSS 
elements. 

In Fig. (top panel) we plott the WmstQ-mst)^ for 
the entire sample of 79,183 galaxies situated at distances 
D < 420/i~ 1 Mpc where, as is seen from Fig. [2 the impact 
of the selection effect is still moderate. The error bars show 
the scatter of measurements for the four subsamples. For 
each sample, Ni, N 2 , Si, S 2 , we have 

(Imst) = 2.5, 2.6, 2.3, 2.0 h^Mpc . (6) 

These variations demonstrate the differences in the sample 
properties (cosmic variance). 

Notice in Fig. El that the Wmst(Imst) is well fit by 
a superposition of Rayleigh (at Imst < (Imst)) and ex- 
ponential (at Imst > (Imst)) functions. This confirms 
results discussed in Doroshkevich et al. (2000; 2001) in 
respect to the high degree of galaxy concentration within 
the population of high density rich wall-like structures 
and less rich filaments. However, as was noted in the same 
papers, with this approach the approximate separation of 
wall-like and filamentary structure elements can be per- 
formed only statistically. This is because the high density 
part of the PDF described by the Rayleigh function in- 
cludes high density clouds situated in both filaments and 
walls. The exponential part of the PDF is related mainly 
to the filamentary component. 

3.3. High and low density regions 

The methods for an approximate statistical decomposition 
of a sample into subsamples of wall-like structures and 



Whdr = (1 ± 0.18) W R , Imst < 1M{I MS t) , (7) 

where Wr is the Rayleigh function JSJ. Larger difference 
between observed and expected PDFs for larger Imst in- 
dicates that the selected sample of HDRs includes some 
fraction of objects, ~ 10%, which can be related to the 
filamentary component with the exponential PDF. 

For the LDRs, the WmStQ) is plotted in Fig. (bot- 
tom panel). For small edge lengths, I < (Imst), it also 
fits well to a Rayleigh function indicating that ~ 60% of 
LDR galaxies are concentrated within less massive 3D (el- 
liptical) and 2D (sheet-like) clouds. This result confirms 
the strong disruption of filaments to a system of high den- 
sity clouds. For larger edge lengths, however, the LDR 
Wmst(1) appears to be closer to an exponential function, 
indicating that according to this criterion the spatial dis- 
tribution of the remaining ~ 40% of LDR galaxies is simi- 
lar to a ID Poissonian one which is typical for filamentary 
structures. 

The mean edge lengths, (Imst), found for HDRs and 
LDRs in samples Nx, N 2 , Si, S 2 are 

(Imst) = 1.2, 1.3, 1.2, 1.3 h~ l Upc , (8) 

(Imst) = 2.8, 2.7, 2.6, 2.8 h^Mpc . (9) 

These values differ by about a factor of two from each 
other, indicating that, as is seen from the difference 
in the mean density within HDRs and LDRs elements is 
about an order of magnitude. Of course, the volume aver- 
aged density of LDRs is still less. 
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Fig. 6. Mass functions of structure elements, / m (e), e = 
Ltr/L S um for the structure elements selected within HDRs 
(top panel, solid and dashed lines) and within LDRs (bot- 
tom panel, solid and dashed lines). 

3.4. Morphology of the structure elements 

Within so defined HDRs and LDRs themselves we can 
extract with the MST technique subsamples of structure 
elements for various threshold overdensities. We can then 
suitably characterize the morphology of each structure el- 
ement by comparing the sum all edge lengths within its 
full tree, L sum , with the sum of all edge lengths within 
the tree's trunk, L tr , which is the longest path that can 
be traced along the tree without re-tracing any steps. The 
ratio of these lengths 

e = Ltr I Lsum. (10) 

suitably characterizes the morphology of the LSS ele- 
ments. 

For filaments, we can expect that the lengths of the 
full tree and of the trunk are similar to each other, e *~ 1, 
whereas for clouds and walls these lengths are certainly 
very different and e < 1. This approach takes into ac- 
count the internal structure of each element rather than 
the shape of the isodensity contour bounding it, and in this 
respect it is co mplementary to the Minkowski Functiona l 
technique fe.g. lSchmalzing et al.lll999llSheth et aljEoO^) . 

However, even this method cannot discriminate be- 
tween the wall-like and 3D (elliptical) clouds and those 
rich filaments having many long branches for which again 
e < 1. This means that both the PDF of this ratio, W(e), 
and the corresponding mass function, / m (e), are contin- 
uous functions and the morphology of structure elements 
can be more suitably characterized by the degree of fil- 
amentarity and 'wall-ness'. This also means that we can 
only hope to distinguish statistical differences between the 



morphologies of structure elements in HDRs and the mor- 
phologies of structure elements in LDRs. 

The selection of clusters within HDRs and LDRs 
was performed for two threshold linking lengths, ri n k = 
2. & 2.4/i~ 1 Mpc for HDRs, and r lnk = 3.2 & SM^Mpc 
for LDRs. These values are larger than the mean edge 
lengths and characterize the LSS elements with interme- 
diate richness when the measured difference between the 
walls and filaments is maximal. As was noted above, for 
lower linking lengths this method characterizes mainly the 
internal structure of the LSS elements while for large link- 
ing lengths filaments percolate and form the joint network 
with again t < 1. 

The distribution functions of the ratio, W(e), are found 
to be close to Gaussian with (e) rs 0.5 & 0.70 for HDRs 
and LDRs, respectively. The mass functions, f m (e), plot- 
ted in Fig.[fJ]for the same linking lengths are shifted to the 
left (for HDRs) and to the right (for LDRs) with respect 
to the middle point. Differences between these functions 
found for smaller and larger linking lengths illustrate the 
impact of the percolation process and disruption of the 
LSS elements. 

These results verify the objective nature of the differ- 
ences in the structure morphologies in HDRs and LDRs. 

4. Statistical characteristics of filaments 

Theoretical characteristics of the LSS elements relate to 
the dominant dark matter (DM) component while the ob- 
served galaxy distribution relates to the luminous matter 
which represent only ~ 3 - 5% of the mean density of 
the Universe. Spatial distribution of dark and luminous 
matter is strongly biased. None the less, some observed 
characteristics of filamentary components of the LSS can 
be compared with both available theoretical expectations 
and characteristics obtained for simulated DM distribu- 
tions. 

The most interesting ones are the PDF of the linear 
density of filaments measured as the mass or number of 
objects per unit length of filament, ~Sfu- The other is the 
mean surface density of filaments (Xfu, defined as the mean 
number of filaments intersecting a unit area of arbitrary 
orientation. Both characteristics depend upon the thresh- 
old linking length, r; n fc, used for the filament selection 
and upon the threshold richness of filaments. However, 
both characteristics are independent from the small scale 
clustering of matter within filaments. 

Comparison of these characteristics of filaments for ob- 
served and simulated catalogues allows one to test the cos- 
mological model used. However, the connection of quan- 
titative characteristics of filaments with the initial power 
spectrum is complex and these measurements cannot yet 
be used for estimates of the power spectrum. 

4.1. Linear density of filaments 

As was found in DD99 and DD02, for the CDM-like initial 
power spectrum the PDF of the filament linear density can 
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are plotted by solid lines. 



by solid lines. 



be written as follows: 

N fa dXfu « exp(-J3E /il /(E /i ,)) cE fa . (11) 

However, poorer filaments cannot actually be selected in 
either simulated or observed catalogues. Hence, even for 
high resolution simulations the measured PDF is well fit- 
ted by the relation 

Nfu « a erf 4 [ai(a; - x )] exp[-y/a 2 (x - x )} (12) 

where x = T, fu / (E, fu) , Xq ~ 0.35, ai ~ 2 — 2.5 and a 2 ~ 
30-40 (DD02) and (£/«) - 3/r /nfc . Here the cutoff of the 
PDF at x — xq reflects the limited resolution of simulated 
matter distribution. 

These results can be compared with ones obtained for 
for the DR1 and the mock catalogue l|Cole et al.l Il998^ 
which allows one to estimate the impact of the selection 
effect for the measured linear density of filaments. In both 
the cases, the linear density of objects was measured by 
the ratio of the number of points and the length of the 
MST for each filaments of the sample. 

For the mock catalogue the PDFs of the linear density 
of filaments are plotted in Fig. \7\ for three linking length 
(with Sthr ~ 1, 0.7 & 0.5). These PDFs are well fitted by 
expression H12fl for parameters xq « 0.6, a% w 2.5, a 2 ~ 
80 — 90. For all linking lengths we get (£ fu)ri n k — 2.25 ± 
0.04. For the observed DR1 catalogue the same PDFs are 
plotted in Fig. [S] also for three linking length (with Sthr ~ 
50, 5 & 0.55). They are well fitted by the same expression 
l|12l) for parameters xo sa 0.5, ai sa 2.5, a 2 ~ 70 — 80. For 
this sample we get (E/ij)n n fe = 2.2 ± 01. 

These results show that, in all the cases, the mea- 
sured PDFs are well fit by the same expression (fL^l 



which coincides with the theoretically expected one Ijllfl 
at T,fu > (Efu). The mean linear density, (£/iz), clearly 
depends upon the linking length used for filament selec- 
tion. The selection effect increases the product (£/jz)j"/ n fc 
by ~ 1.5 times as compared with results obtained for the 
DM simulations. 

These results show that the observed galaxy distribu- 
tion nicely represents the expected and simulated ones. 
The results also indicate that the general properties of 
filaments are consistent with a CDM-like initial power 
spectrum. 

4.2. Typical size of the filamentary network 

Due to the complex shape of the network of filaments 
spanning the LDRs any definition of the typical size of 
a network cell is m e arly c onvenient. As was discussed in 
iDoroshkevich et "all l|200lh , two definitions seem to be the 
most objective. One is the mean free path between fila- 
ments along a random straight line. The other is the mean 
distance between branch points of the tree along the trunk 
of selected filaments. The second definition tends to yield 
cell sizes that are typically a factor of 1.5 smaller than 
those yielded by the mean-free-path definition. 

Theoretical estimates of this size are uncertain be- 
cause it strongly depends upon the sample of selected 
filaments (DD02). This means that this characteristic 
strongly depends upon the catalogue used. Moreover, fil- 
aments are connected to the network only for larger link- 
ing lengths; thus the typical measured cell size depends 
also upon the threshold linking length used. Hence, for 
the LCRS the mean free path between filaments with a 
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Fig. 9. Distribution functions, N, for the distance between 
branch points along a trunk for filaments selected in LDRs 
(thick solid and dashed lines). Fit (|T3)) is plotted by thin 
solid line. 



variety of richness was estimated in iDoroshkevich et "all 
(2001) as ~ 13 — 30ft, 1 Mpc. The mean distance between 
branch points of the tree along the trunk was estimated 
as pa 10ft. _1 Mpc and it rapidly increases with the linking 
length used owing the progressive percolation of filaments 
and formation of the joint LSS network. 

Here with a richer sample of filaments we can also es- 
timate the PDFs of the cell sizes measured by the dis- 
tance between branch points of the tree along the trunk. 
These PDFs, N(lb r ), are plotted in Fig.OHfor two linking 
lengths, ri n k = 1-8 & 3.6/i _1 Mpc, which correspond to the 
threshold overdensities Sthr = 0.66 & 0.5. These PDFs are 
roughly fitted by expression 

N(l br ) w 270x 4 - 5 exp(-9.1a;), x = l br /{l br ) , (13) 
The measured mean distance between branch points, 
(l br ) « 4.7 & 11.9/i -1 Mpc. 

are close to those obtained in IDoroshkevich et all l)200lh 
and IDoroshkevich et alJ 



1996) and those cited above. For 
smaller n n k this estimate is decreased because of the dom- 
ination of short filaments which are not yet connected to 
the network. 



5. Parameters of the wall-like structure elements 

The statistical characteristics of obse rved walls were first 
measured using the LCRS and DURS <)Doroshkevich et alJ 
l2000U200ll) . The rich sample of walls extracted from the 
SDSS DR1, however, permits more refined estimates of 
these characteristics. As was discussed in Sec. 3.2 walls 



dominate the HDRs, and thus these subsamples of galaxies 
can be used to estimate the wall properties. 

The expected characteristics of wa lls and methods 
of the ir measurement were discussed in iDemianski et all 
(2000) so here we will only briefly reproduce the main 
definitions. It is important that these characteristics can 
be measured independently in radial and transverse direc- 
tions, which reveals the strong influence of the velocity 
dispersion on other wall characteristics. 

5.1. Main wall characteristics 

Main characteristic of walls is their mean dimensionless 
surface density, (q w ), measured by the number of galaxies 
per Mpc 2 and normalized by the mean density of galaxies 
multiplied by a coherent length of the initial velocity field 
(DD99; DD02) 

l v re 33/i, _1 Mpc (0.2/r), T = Q m h , (14) 

where Q m is the mean matter density of the Universe. For 
Gaussian initial perturbations, the expected probability 
distribution function (PDF) of the surface density is 



1 



1 



N th {q w ) 
(q th ) = 8(0.5 + I/ttJt* 



exp 



q w 
'8t1 



erf 



it I • < 15 > 



6.55r„ 



This relation links the mean surface density of walls with 
the dimensionless amplitude of perturbations, r m , 



= \/ (q w )/6.55. 



(16) 



which can be compared with those measured by other 
methods (DD02, Sec. 8.2). 

Other important characteristics of walls are the mean 
velocity dispersion of galaxies within walls, (w w ), the 
mean separation between walls, (D aep ), the mean over- 
density, (S), and the mean thickness of walls, (h). The 
mean velocity dispersion of galaxies, (w w ), can be mea- 
sured in radial direction only whereas other wall character- 
istics can be measured both radially and along transverse 
arcs. Comparison of the wall thickness and the overden- 
sity, (h) and (S), measured in transverse (t) and radial (r) 
directions, illustrates the influence of the velocity disper- 
sion of galaxies on the observed wall thickness. 

The velocity dispersion of galaxies within a wall w w 
can be r elated to the radial thi ckness of the wall by this 
relation ijDemianski et al JI20 001 : 



h r = Vl2H 1 w v 



(17) 



For a relaxed, gravitationally confined wall, the mea- 
sured wall overdensity, surface density, and the velocity 
dispersion are linked by the condition of static equilib- 
rium. Consider a wall as a slab in static equilibrium, and 
this slab has a nonhomogeneous matter distribution across 
it. We can then write the condition of static equilibrium 
as follows: 

2 ttG/z 2 



(p)s 



a, 



8 6 



■{Hol v q w ) 0$ . 



(18) 
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Here \i = (p)l v q w is the mass surface density of the wall 
and the factor 0$ ~ 1 describes the nonhomogeneity of 
the matter distribution across the slab. Unfortunately, for 
these estimates we can only use the velocity dispersion and 
overdensity measured for radial and transverse directions. 
Hence, the final result cannot be averaged over the samples 
of walls. 

5.2. Measurement of the wall characteristics 

The characteristics of the walls can be measured 
with the two paramet er core-sampling approach 
l|Doroshkevich et alJ [l996) applied to the subsample 
of galaxies selected within HDRs. With this method, all 
galaxies of the sample are distributed within a set of 
radial cores with a given angular size, 9 C , or within a set 
of cylindrical cores oriented along arcs of right ascension 
with a size d c . All galaxies are projected on the core axis 
and collected to a set of one-dimensional clusters with a 
linking length, lu n k- The one-dimensional clusters with 
richnesses greater than some threshold richness, N m i n , 
are then used as the required sample of walls within a 
sampling core. 

Both the random intersection of core and walls and the 
nonhomogeneous galaxy distribution within walls lead to 
significant random scatter of measured wall characteris- 
tics. The influence of these factors cannot be eliminated, 
but it can be minimized for an optimal range of param- 
eters 9 Cl d Cl kink and N m i n . Results discussed below are 
averaged over the optimal range of these parameters. 

For the measurement of wall characteristics in the ra- 
dial direction four samples of HDRs galaxies were used 
in each field of the DR1 catalogue. One of these sam- 
ples was selected as was discussed in Sec. 3.3, three other 
samples were selected from the catalogues already cor- 
rected for radial selection effects (Sec. 2.1) with the same 
threshold overdensity Sthr — 1 and for HDRs contain- 
ing ~ 43%, 50%, and 56% of all galaxies in the field. In 
all the cases, the wall parameters were measured in real 
space for the selected samples of the HDRs. The mean 
wall properties were averaged over four radial core sizes 
(9 C = 2°, 2.25°, 2.5° and 2.75°) and for six core-sampling 
linking lengths (2/i _1 Mpc< l lin k < 4.5/i" 1 Mpc). Final av- 
eraging was performed over all sixteen samples and over 
all 9 C and l Unk . 

Due to the complex shape of the fields S2 and N2, the 
measurements of wall characteristics in the transverse di- 
rection were performed for the fields S% and Ni only. The 
mean wall properties were averaged over four core diam- 
eters (d c — 6.0, 6.5, 7.0, and 7.5 h~ 1 Mpc) and five core- 
sampling linking lengths (2/i _1 Mpc < ln n f- < 4./i~ 1 Mpc). 

5.3. Measured characteristics of walls 

The mean radial and tranverse wall properties for all fields 
are listed separately in Table 1. Characteristics obtained 
by averaging over all samples are compared with those 



from the DURS and LCRS and wit h those from mo ck cat- 
alogues simulating the SDSS EDR l|Cole et al.ll998f) . Both 
DM simulation and mock catalogues are prepared for the 
ACDM cosmological model with Q m = 0.3, fl\ — 0.7 and 
with the amplitude of perturbations ffs = 1.05 that excee d 
the now accepted value er s = 0.9±0.1 IjSpergel et al . 2003). 
This excess of the amplitude increases both the measured 
r m and (q) by ~ 10% and ~ 20%, respectively. 

The richness and geometry of these catalogues are 
strongly different. Thus, DURS is an actual 3D cata- 
logue but it contains ~ 2 500 galaxies at the distance 
D < 250/i -1 Mpc and its representativity is strongly lim- 
ited. The LCRS include ~ 21 000 galaxies at the distance 
D < 450/i _1 Mpc but they are distributed within six thin 
slices that again distorts the measured wall characteris- 
tics. Moreover, both catalogues include a small number of 
walls with large scatter in richness. This leads to a signif- 
icant scatter of measured characteristics of walls for these 
catalogues which is a manifestation of well known cosmic 
variance. Only with an actually representative catalogue 
such as the the SDSS can this effect can be suppressed. 

The difference between the mean wall surface densities 
measured for ~ 15-20% of samples reflects real variations 
in wall properties for a limited portion of the samples. 
However, the scatters of mean values listed in Table 1 
partially include the dispersions depending on the shape 
of their distribution functions. The actual scatter of the 
mean characteristics of walls averaged over all samples 
listed in Table 1 is also < 10-12%. 

The amplitude of initial perturbations characterized 
by values r m for the richer sample measured in the radial 
direction is 

rw (0.53 ± 0.05)VT = (0.24 ± 0.02) (19) 

This is quite consistent with estimates found for simula- 
tions. Differences between this value and T m measured in 
the transverse direction for the DR1, the LCRS and the 
DURS demonstrate the impact of the representativity of 
the catalogue used. 

The difference between the wall thickness measured in 
the radial and transverse directions, h r and h tl indicates 
that, along a short axis, the walls are gravitationally con- 
fined stationary objects. Just as with the 'Finger of God' 
effect for clusters of galaxies, this difference characterizes 
the gravitational potential of compressed DM rather than 
the actual wall thickness. The same effect is seen as a dif- 
ference between the wall overdensities measured in radial 
and transverse directions. 

The difference between the wall thicknesses is com- 
pared with the velocity dispersions of galaxies within the 
walls, (w w ). Clusters of galaxies with large velocity dis- 
persions incorporated in walls also increase the measured 
velocity dispersion. The correlation between the wall sur- 
face density and the velocity dispersion confirms the relax- 
ation of matter within walls. This relaxation is probably 
accelerated due to strong small scale clustering of matter 
within walls. 
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Table 1. Wall properties in observed and simulated catalogues 



sample 


Ngal 


<?»}/r 




{St) 


<M 

h^Mpc 


<M 


(w w ) 












fc -1 Mpc 


km/s 


/i _1 Mpc 










radial cores 












Nl 


41217 


2.01 ±0.21 


0.55 ±0.03 1.3 




10.2 ± 1.9 


- 


295 ± 55 


69 ± 14 


N2 


25 935 


1.61 ±0.16 


0.49 ±0.03 1.1 




10.5 ±2.0 


- 


303 ± 57 


78 ± 17 


SI 


13 215 


2.21 ±0.19 


0.58 ±0.03 1.4 




10.5 ± 1.7 


- 


303 ± 50 


67 ± 14 


S2 


12 585 


1.79 ±0.33 


0.51 ±0.04 1.3 




9.0 ± 1.8 


- 


260 ± 50 


83 ±20 






transverse cores for the SDSS DR1 










Nl 


16 883 


2.47 ±0.51 


0.61 ±0.07 


3.9 




4.3 ±0.8 


- 


65 ± 11 


SI 


13 215 


2.29 ±0.64 


0.58 ± 0.08 


3.9 




4.0 ±0.8 


- 


58 ± 11 








observed samples 










SDSS (radial) 


92 952 


1.91 ±0.32 


0.53 ±0.05 1.3 




10.1 ± 2.0 


- 


291 ± 56 


74 ± 17 


SDSS (transverse) 


29311 


2.42 ±0.67 


0.60 ±0.08 


3.5 




4.9 ± 1.3 


- 


64 ± 14 


LCRS 


16 756 


2.51 ± 0.9 


0.62 ±0.10 3.0 


7.4 


8.6 ±0.8 


2.8 ±0.7 


247 ± 48 


60 ± 10 


DURS 


2 500 


2.23 ±0.6 


0.58 ±0.08 1.7 


6.5 


9.7 ± 1.8 


4.9 ± 1.2 


280 ± 52 


44 ± 10 




mock catalogues in real and redshift spaces 


for the model with T = 


3.2 






redshift 


98 828 


2.7 ±0.5 


0.63 ±0.06 1.8 


3.8 


11.8 ±2.1 


6.5 ± 1.4 


338 ± 65 


50 ± 10 


real 


98 828 


2.1 ±0.4 


0.57 ±0.06 4.3 


4.6 


4.8 ± 1.0 


4.2 ± 1.0 


305 ± 47 


50 ± 10 




DM catalogue in 


real space for the model with F 


= 0.2 








real 


7.1 • 10 6 


2.5 ±0.4 


0.63 ±0.04 2.7 




4.9 ±0.5 




245 ± 30 


52 ±5 



o 




Fig. 10. The PDFs of dimensionless surface density of 
walls, N m (q/ (<?)), for walls selected in the DM simulation 
(top panel) and the mock catalogue in real (middle panel) 
and redshift (bottom panel) spaces. Fits (|15fl are plotted 
by solid lines. 



o 
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Fig. 11. The PDFs of observed dimensionless surface den- 
sity of walls, N m (q/ (q)), for walls selected in four regions 
of the SDSS galaxy catalogues. Theoretically expected fits 
are plotted by solid lines. 



The measured PDFs of the surface density of walls are 
plotted in Figs. EH and ^2 f° r the simulated DM distri- 
bution, mock catalogues in real and redshift spaces, and 
four observed samples of the SDSS DR1. These are nicely 
fit by the expected expression fT5|l . Thus, for simulated 



samples we have, respectively, 

N m = (1 ± 0.1)N th , (q)/(q th ) = 0.78, 0.87, 0.87 , (20) 

where Nth an d (qth) are given by ijTo)) . For the two south- 
ern samples of the SDSS DR1 we get: 

N m = (1.1 ± Q.2)N th , (q)/(q th ) = 1.07 , (21) 
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Fig. 12. The PDFs, N w (w/(w)), of reduced velocity dis- 
persion within wall 124|) for walls selected in four regions 
of the SDSS galaxy catalogues. Gaussian fits are plotted 
by solid lines. 

and for the two northern samples of the SDSS DM we 
get: 



N m = (1. ± 0.1)N th , (q)/(q th ) =0.97, 



(22) 



These results verify that, the observed walls represent re- 
cently formed Zel'dovich' pancakes. 

Using the measured mean wall overdensity in the 
transverse direction listed in Table 1 we have for the pa- 
rameter 0$ introduced in equation E 



(S) 0.3 

T~ o 



1.1. 



(23) 



which is also consistent with the expected values for re- 
laxed and stationary walls. 

As was proposed in iDemiahski et alJ l)2000h we can 
discriminate between systematic variations in the mea- 
sured velocity dispersion due to regular variations in the 
surface density along the walls (Fig. Illfl and the random 
variations in the velocity dispersion which integrates the 
evolutionary history of each wall. Indeed, along a shorter 
axis, for gravitationally bound and relaxed walls we can 
expect that 



oc q 2 /S cx <5 7 



w w oc q 



1-1/7 



Here we assume that the distribution of DM component 
and galaxies can be approximately described by the poly- 
trop ic equation of state wit h the power index 7 ss 5/3 
- 2. iDemianski et alJ ((2003) suggest for consideration a 
reduced velocity dispersion, u) w , 
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Fig. 13. The PDFs, N sep (D sep / (D sep )), of observed wall 
separations for walls selected in four regions of the 
SDSS galaxy catalogues. Theoretically expected for the 
Possonian distribution exponential fits are plotted by solid 
lines. 



corrected for variations of the wall thickness. For this func- 
tion the systematic variations of w w are essentially sup- 
pressed and, in most respects, it is similar to the entropy 
of compressed matter. It integrates the action of random 
factors in the course of wall formation. Hence, for this 
function the Gaussian PDF, N^, can be expected. Indeed, 
this PDF plotted in Fig. HI for four samples of the SDSS 
DR1 is quite similar to a Gaussian function with a stan- 
dard deviation <j u ps 2.25(w). These results show a large 
scatter of evolutionary histories of observed walls. 

Note that, for all the samples listed in Table 1, the 
mean wall separation, (D sep ), is close to twice that of the 
coherent length of the initial velocity field, 



(D 



sep/ 



2 lq; 



(25) 



oj w = I ln(w w q w Pw )\, p w ps 1 - I/7 fa 0.5. 



(24) 



for the low density cosmological models with T ps 0.2 l|14fl. 
These results coincide with the estimates of the matter 
fraction, ~ 50%, accumulated within walls. Due to the 
large separation of walls, the correlations of their posi- 
tions is small and a random ID Poissonian PDF of the 
separation can be expected. These PDFs are plotted in 
Fig. 1131 together with the exponential fits. 

Finally, we would like to draw attention to the fact that 
all measured properties of these walls are quite consistent 
with a CDM-like initial power spectrum and Gaussian 
distribution of perturbations. 
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6. Possible rich clusters of galaxies 

The SDSS DM also contains a number of galaxy com- 
plexes of various richnesses which can be extracted by 
means of the MST technique. Due to the large velocity dis- 
persion of galaxies within clusters and the strong 'Finger 
of God' effect, this extraction must be performed using 
different threshold linking lengths in the radial (l r ) and in 
the transverse (l t ) directions. This is not unlike how group 
catalogs are extracted from redshift s urveys using con- 
ventional 'friends-of- frien ds' algorithms ijHuchra fc Gellcr 
Il982t iTucker et al.ll2000l) . 

We performed this cluster-finding in two major steps. 
First, we projected the observed samples onto a sphere 
of radius R = lOO/i^Mpc with a random scatter 
±0.5/i -1 Mpc and extracted a set of candidate clus- 
ters from this catalog using a linking length of r t = 
0.3/i^ 1 Mpc. Second, we applied a radial linking length 
of r> = 3/i _1 Mpc to these candidate clusters using their 
real 3D coordinates. In this second step, we also employed 
the threshold richness, N mem = 10, for our final samples 
of possible rich complexes. Having extracted these proba- 
ble rich "clusters", we calculated a distance-independent 
measure of their richnesses by correcting their observed 
richnesses N mem for radial selection effects using equation 
(1); we call this corrected richness N se i. Further discrimi- 
nation of the "clusters" can be performed using their size 
in transverse and radial directions as well as other char- 
acteristics of selected galaxies. 

For the threshold parameters used 20, 12, 43 
and 10 possible rich clusters were selected from the 
Si, S2, Ni & N 2 samples, respectively. The majority of 
these "clusters" are embedded within richer walls. Let 
us remember that these are possible rich clusters of 
galaxies and to confirm that they are physical poten- 
tial wells, it is necessary to check for diffuse x-ray emis- 
sion. However, compar ison with the list of NORAS survey 
ijBohringer et al.ll2000|) shows that only ~ 32 of them can 
be roughly related to the x-ray sources. This means that 
the method of cluster identification must be essentially 
improved. 

7. Mass function of the structure elements 

The richness of the SDSS DR1 allows one to extract sev- 
eral different sets of high density clouds and structure el- 
ements with various overdensities within the HDRs and 
LDRs and to find their mass function. These results can 
be directly compared with the theoretical expectations of 
DD02. However, the richness of the SDSS DR1 does not 
yet allow one to estimate quantitatively the divergence be- 
tween the expected and observed mass functions plotted 
in Fig. H3 

Two samples of high density galaxy groups and two 
samples of unrelaxed structure elements - walls and fila- 
ments - were selected separately within HDRs and LDRs 
introduced in Sec. 3.3 for a threshold richness N mem > 5. 
Since the velocity dispersions in groups are expected to be 



much smaller than those in rich clusters, we select these 
samples of structure elements using the simpler method 
described in Sec. 4 rather than the two-step approach de- 
scribed in Sec. 6. The richness of each cluster was corrected 
for radial selection effects using the selection function in- 
troduced in Sec. 2.1. 

The main properties of the selected clouds are listed in 
Table 2, where ri n k and (6) are, respectively, the thresh- 
old linking length and the mean overdensity of clouds, 
fg a i is the fraction of galaxies from the total (combined 
HDR+LDR) sample of galaxies within the selected clouds, 
N c i is the number of clouds while (M) and (M se i) are the 
observed and corrected for the selection effect mean rich- 
ness of individual clouds. 

Results listed in Table 2 illustrate the influence of 
environment on the properties of high density clouds. 
In particular, in spite of the approximately equal num- 
ber of galaxies in HDRs and LDRs, ~85% and ^70% 
of the high density clouds selected with linking lengths 
fink = 0.8 & 1.2/i _1 Mpc are situated within the HDRs 
and accumulate ~90% and ~80% of galaxies. At linking 
length rink — 1.8/i _1 Mpc the numbers of clouds selected 
within HDRs and LDRs are comparable but again ^80% 
of galaxies related to these clouds are concentrated within 
HDRs. At the largest linking lengths listed in Table 2 es- 
sential differences are seen only for the mean masses and 
number of structure elements selected within HDRs and 
LDRs. These differences are enhanced by the influence of 
the selection effect which is stronger for the LDRs. 

The mass functions for these samples are plotted in 
Fig. As was shown in DD02, in Zel'dovich theory and 
for the WDM initial power spectrum the dark matter 
mass function of structure elements is independent of their 
shapes and, at small redshifts, it can be approximated by 
the functions 

xN(x)dx = 12.5k Z ax 2/3 exp(-3.7x 1/3 ) erf (x 2/3 )dx , (26) 
xN(x)dx = 8.k za x 1/2 exp(-3.1a; 1/3 ) erf {x 3/i )dx . (27) 
M 

The expression <|26|) relates to clouds which have become 
essentially relaxed and static by z = 0, and the expres- 
sion (|27|l relates to richer, unrelaxed filaments and walls 
which are still in the process of collapse. Here, kza ~ 
1.5-4 and /iza ~ 0-8 - 1.3 are fit parameters which 
take into account the incompleteness of selected samples of 
clouds for small and large richnesses; this incompleteness 
changes both the amplitude and mean mass of the mea- 
sured clouds. Comparison with simulations (DD02) has 
shown that these relations fit reasonably well to the mass 
distribution of DM structure elements. 

These relations are similar to the mass function from 
the Press- Schechter formalism, 

xN PS (x)dx = -^L C 1/6 exp(-^ 3 )dx , (28) 
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Fig. 14. Mass functions of galaxy clouds, N m ■ N mem / (N mem )) , selected in HDRs (left panels) and 
LDRs (right panels) for four threshold linking lengths. Theoretical fit (|26[1 for relaxed clouds (ri n k = 
0.8&1.2/i _1 Mpc) and fit l|27(l for unrelaxed clouds (rink > 1.8ft. -1 Mpc) are plotted by solid lines. 



Table 2. Parameters of groups of galaxies selected in 
HDRs and LDRs after correction for the selection effect. 



r lnk h i Mpc 


(S) 


j gal 


N cl 


(M) 


(M sel ) 






HDR 








0.8 


140 


0.1 


1193 


6.8 


12 


1.2 


40 


0.24 


1961 


9.5 


25 


1.8 


10 


0.39 


1731 


17.7 


61 


2.6 


3 


0.47 


614 


61.5 


228 






LDR 








0.8 


220 


0.014 


220 


4.8 


14 


1.2 


61 


0.05 


817 


5.3 


18 


1.8 


16 


0.15 


1805 


6.4 


36 


3.6 


3 


0.36 


2 295 


12.5 


106 



£ = 1.785/vz = 1.785/i ps M/(M) , 

despite the fact that they use different assumptions about 
the process of cloud formation and the shape of the formed 
clouds. Here again the fitting parameters Hps and [ips 
take into account the incompleteness of measured sample. 
However, this expression relates to the CDM-like power 
spectrum without small scale cutoff which is linked, for 
example, with the finite mass of DM particles, and without 
correction for the survival probability. So, it describes only 
the massive part of the mass function. 

Relations (I26[) . & (|27|) characterize the mass distri- 
bution of dark matter clouds associated with the ob- 
served galaxy groups and massive structure elements. 
They are closely linked with the initial power spectrum 
and fit reasonably well the observed mass distribution. 
For N se i < (N se i) the incompleteness of the sample of se- 



lected clouds leads to rapid drops in the observed mass 
functions as compared with theoretical expectations. For 
filaments selected within LDRs at r;„fc < 2 the deficit of 
richer clouds is enhansed by the method of filament sep- 
aration. However, for the largest linking lengths, r;„fc = 
2.6 & 3.6/i _1 Mpc, where the incompleteness and other 
distortions are minimal, the observed mass distribution is 
quite consistent with theoretical expectations. 

8. Summary and discussion 

Statistical analysis of large galaxy redshift surveys allows 
us to obtain the quantitative characteristics of the large 
scale galaxy distribution, which in turn can be related to 
the fundamental characteristics of the Universe and the 
processes of structure formation. The large homogeneous 
data set compiled in the SDSS DR1 also permits us to 
checking the results from analysis of the LCRS and the 
DURS and to obtain more accurate and more represen- 
tative estimates of the main basic characteristics of the 
Universe. 

The spatial galaxy distribution for the Si samples is 
plotted in Fig. 15; galaxies in HDRs are highlighted. 

8.1. Main results 

The main results of our investigation can be summarized 
as follows: 

1. The analysis performed in Sec. 3 with the MST tech- 
nique confirms that about half of galaxies are situated 
within rich wall-like structures and the majority of the 
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demonstrate the impact of environment and interac- 
tion of small and large scale perturbations. 
5. We found the typical cell size of the filamentary net- 
to be ~ 10ft. -1 Mpc (Sec. 4). This estimate is con- 
sent with the one obtained pr eviously for the LCRS 
j Doroshkevich et al.l[l996tl200l|) . 
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Fig. 15. Spatial distribution of all galaxies in the Si sam- 
ple (points) and within the HDRs (thick points). 



remaining galaxies are concentrated within filaments. 
This result confirms that the filaments and walls are 
the main structure elements in the observed galaxy 
distribution. Quantitative characteristics of walls and 
filaments presented in Sections 3, 4 & 5 validate this 
division of the LSS into these two subpopulations. 

2. The main characteristics of wall-like structure ele- 
ments, such as the overdensity, separation distance be- 
tween walls, wall thickness, and the velocity dispersion 
within walls, were measured separately for radial and, 
for the SDSS DR1 equatorial stripes, in transverse di- 
rections. Comparison of these characteristics demon- 
strates that the walls are approximately in static equi- 
librium, that they are relaxed along their shorter axis, 
and that their observed thickness in the radial direc- 
tion is defined by the velocity dispersion of galaxies. 

3. The PDF of the wall surface density is consistent 
with the simulated one and with that predicted by 
Zel'dovich theory for Gaussian initial perturbations 
and CDM-like initial power spectrum. The measured 
amplitude of perturbations coincides with that ex- 
pected for the spatially flat ACDM cosmological model 
with ft a ~ 0.7 and O m « 0.3. These results demon- 
strate that the spatial galaxy distribution traces the 
dark matter distribution nicely. 

4. The mass distributions of groups of galaxies, filaments 
and walls selected with various threshold overdensi- 
ties nicely fit the joint mass function consistent with 
the expectations of Zel'dovich theory. Characteristics 
of these groups of galaxies listed in Table 2 clearly 



8.2. Characteristics of the LSS 

Results obtained in Sec. 3 with the MST technique demon- 
strate that at least 80 - 90% of galaxies in the SDSS 
DR1 are concentrated within the filaments and walls with 
a var ety of richness and overdensity forming the LSS. 
Comparison of the mean edge lengths (OIH1 and ^ for the 
full sample, HDRs and LDRs with the mean density of the 
sample © shows that the LSS occupies only roughly half 
of the volume. It also shows that the mean overdensity of 
the walls is larger than 10. The same results show that the 
filaments are formed by a system of high density groups 
' Of ' gal axies which contain ss 30% of all galaxies. For only 
« 25%3of galaxies do we see the actual one dimensional 
PDF of the MST edge lengths (Fig.EJ. 

There is continuous distribution in morphology and 
richness of the LSS elements. It can be roughly described 
as a system of walls randomly distributed in space with 
mean separations of « (50 ± 10)/i _1 Mpc and a filamen- 
tary network connected to the walls. The typical size of 
lower density regions, or voids, situated between filaments 
and walls is estimated in Sec. 4 as approximately 10-12 

These results agree with those obtained for the mock 
catalogues simulating the galaxy distribution and indicate 
that galaxies nicely trace the spatial distribution of dom- 
inated dark matter. 



8.3. Properties of walls and parameters of the initial 
power spectrum 

Walls and filaments are the largest structure elements ob- 
served in the Universe. In contrast to galaxies, their for- 
mation occurs at relatively small redshifts in course of 
the last stage of nonlinear matter clustering and is driven 
by the initial power spectrum of perturbations. Therefore, 
their properties can be successfully desc ribed by the non- 
linear theory of gravitational instability ( ZerDovichll 970'l 
that allows us to link them with the characteristics of the 
initial power spectrum. 

The interpretation of the wa lls as Zel'dovich pan- 
cakes has bee n discussed already in lThomnson & Gregory! 
1)19781) and in lOortl 1)19831) . The comparison of the statisti- 
cal characteristics of the Zel'dovich pancakes for a CDM- 
like initial power spectrum (DD99, DD02) with those for 
observed walls demonstrates that, indeed, this interpre- 
tation is correct and, for a given cosmological model, it 
allows us to obtain independent estimates of the funda- 
mental characteristics of the initial power spectrum. 
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The estimates of the mean wall surface density, (q w ), 
and the amplitude of initial perturbations, (r m ), listed in 
Table 1 are consistent with each other and with those 
found for the LCRS and DURS. They are also close to 
those found for the simul ated DM distrib ution and for 
the mock galaxy catalogs l|Cole et al.lll998F) prepared for 
a spatially flat ACDM cosmological model with Cl\ = 0.7, 
Q m = 0.3 and as = 1.05. As was shown in Sec. 5.3 <(2TJ|) . 

and , the PDFs of both observed and simulated 
wall surface density plotted in Figs. and El coincide 
with those theoretically expected (|15f) for Gaussian initial 
perturbations with a CDM-like power spectrum (DD99; 
DD02). 

Averaging of both (q w ) and (r m ) listed in Table 1 al- 
lows us to estimate the mean values as follows: 

(q w ) = (0.38 ±0.06)(r/0.2), (29) 



r m = (0.24 ± 0.02) A/r/0.2 . (30) 

These values are consistent with the best estimates of the 
same amplitude ijSpergel et al.ll200^l for the ACDM cos- 
mological model with T — 0.2 

erg » 0.9 ±0.1, t« 0.22 ±0.02, (31) 

These results verify that galaxies nicely trace the LSS 
formed mainly by the dark matter distribution and the 
observed walls are recently formed Zel'dovich pancakes. 
They also verify the Gaussian distribution of initial per- 
turbations and coincide with the Harrison - Zel'dovich 
primordial power spectrum. 

Comparison of other wall characteristics measured in 
radial and transverse directions indicate that the walls 
are gravitationally confined and relaxed along the shorter 
axis. The same comparison allows us to find the true wall 
overdensity, wall thickness, and the radial velocity disper- 
sion of galaxies within walls. As is seen from relation l|18[) . 
these values are quite self-consistent. 

8.4. Mass function of structure elements 

The samples of walls, filaments, and groups of galaxies in 
the SDSS DR1 selected using different threshold overden- 
sitics allow us to measure their mass functions, to trace 
their dependence on the threshold overdensity and envi- 
ronment, and to compare them with the expectations of 
Zel'dovich theory. 

This comparison verifies that for lower threshold over- 
densities for both filaments and wall-like structure ele- 
ments, the shape of the observed mass functions is consis- 
tent with the expectations of Zel'dovich theory. However, 
for high density groups of galaxies some deficit of low mass 
groups caused, in particular, by selection effects and en- 
hanced by the restrictions inherent in our procedure for 
group-finding leads to a stronger difference between the 
observed and expected mass functions for M se [ < (M se i). 
The same factors distort the observed mass functions for 
groups selected within LDRs. 



Let us note that mass functions (|26II27I & are 
closely linked with the initial power spectrum. This is 
manifested as a suppression of the PDFs at M se ; < (M se i) 
and is proportional to exp[— (M/(M se i)) 1 ^ 3 ] at M se i > 
(M se i). They differ from the mass function of galaxy clus- 
ters and the probable mass function of observed galaxies 
which are formed on account of multi-step merging of less 
massive clouds and are described by a power law with a 
negative power index at M se i < (M se i) and an exponen- 
tial cutoff oc ex£j-(M/(M seZ )) at M sel > (M sel ) (see e.g. 
ISilk fc Whitelll978l) . 

8.5. Interaction of large and small scale perturbations 

The data listed in Table 2 shows that the majority of 
high density groups of galaxies and the main fraction of 
galaxies related to these groups (up to 80 - 90%) are sit- 
uated within HDRs. These results illustrate the influence 
of environment on galaxy formation and the clustering 
of luminous matter. It also indicates the importance of 
interactions of small and large scale perturbations for the 
formatio n of the observed LSS. T his problem was also dis- 
cussed in lEinasto et all (2003a b). 

These differences arc partly enhanced by the influence 
of selection effects. Indeed, the majority of HDRs are situ- 
ated at moderate distances D ~ 150 — 300/i _1 Mpc where 
this influence is not so strong. On the other hand the LDRs 
also include all galaxies situated in the farthest low den- 
sity regions. However, as was shown in DD02 significant 
differences between the characteristics of clouds separated 
within HDRs and LDRs is seen even for simulated clus- 
tering of dark matter. 

These results are consistent with the high concentra- 
tion of observed galaxies within the filaments and walls 
forming the LSS noted in Sec. 8.2. Bearing in mind that 
the mean density of luminous matter does not exceed 10% 
of the full matter density of the Universe we can con- 
clude that galaxies are formed at high rcdshifts presum- 
ably within compressed regions which are now seen as el- 
ements of LSS. 

A natural explanation for both of these differences and 
the high concentration of galaxies within the LSS elements 
is the interaction of small and large scale perturbations 
when large scale compression of matter accelerates the 
formation of small scale high density clouds. Such interac- 
tions may also explain the existence of large voids similar 
to the Bootes void where formation of galaxies has been 
strongly suppressed. 

8.6. Possible rich clusters of galaxies 

The MST technique generalizes the standard 'friends-of- 
friends' method of the selection of denser clouds and of 
probable clusters of galaxies. However, first attempts of 
such selection presented in Sec. 6 show that there is an es- 
sential difference between the selected high density clouds 
and x-ray sources. The nature of this difference is not yet 
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clear and perhaps it will be eliminated after the intro- 
duction of stronger criteria for the selection of probable 
clusters of galaxies from the survey. 

8.7. Final comments 

The SDSS llYorket alJ 120001: IStoughton et~aT] 120021: 
lAbazaiian et alJ l2003|) and 2dF l|Colless et all l200l[) 
galaxy redshift surveys provide deep and broad vistas with 
which cosmologists may study the galaxy distribution on 
extremely large scales -scales on which the imprint from 
the primordial fluctuation spectrum has not been erased. 

In this paper, we have used the SDSS DR1 to investi- 
gate the galaxy distribution at such large scales. We have 
confirmed our earlier results, based on the LCRS and 
DURS samples, that galaxies are distributed in roughly 
equal numbers between two different environments: fil- 
aments, which dominate low density regions, and walls, 
which dominate high density regions. Although different 
in character, these two environments together form a frag- 
mented joint random network of galaxies - the cosmic web. 

Comparison with theory strongly supports the idea 
that the properties of the observed walls are consis- 
tent with those for ZePdovich pancakes formed from 
a Gaussian spectrum of initial perturbations for a flat 
ACDM Universe (Qa ~ 0.7, O m » 0.3). These results are 
consistent with the e stimate of T = 0.20±0.03 obtained in 
IPercival et all (l200lh for the 2dF Galaxy Redshift Survey 
(see also ISnergel et af]l2003h . 

Such analysis one allows to obtain some important ba- 
sic conclusions regarding the properties and the process of 
formation and evolution of the large scale structure of the 
Universe. With future public releases of the SDSS data 
set, we hope to refine these conclusions. 
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